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Abstract 

We assume cylindrically symmetric stars which begin collapsing by 
dissipating energy in the form of heat flux. We wish to study the effects 

4 

of Carroll-Duvvuri- Trodden- Turner (CDTT) model, f(R) = R + a^, 
on the range of dynamical instability. For this purpose, perturbation 
scheme is applied to all the metric functions, material functions and 
f(R) model to obtain the full set of dynamical equation which control 
the evolution of the physical variables at the surface of a star. It is 
found that instability limit involves adiabatic index T which depends 
on the density profile and immense terms of perturbed CDTT model. 
In addition, model is constrained by some requirement, e.g. positivity 
of physical quantities. We also reduce our results asymptotically as 
/i — > 0, being the GR results in both the Newtonian and post Newto- 
nian regimes. 

1 Introduction 

An increasing attention has been paid to the modification of Einstein-Hilbert 
(EH) action with the lagrangian density y/—gR, where g is the metric tensor 
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and R is the curvature scalar. A simplest modification is to add terms which 
are proportional to ^f^g~R n . It is known that for n > 1, terms help to 
understand the standard cosmology at early times showing de-Sitter behavior 
[1]. For n < 0, these corrections provide a possible gravitational alternative 
to dark energy by self-accelerating vacuum solutions and hence lead to the 
cosmic speed-up [2J. 

The most studied, simplest and definiteness f(R) models in the literature 
are CDTT model f(R) = R + cr^-, super gravity model f(R) = R + aR 2 
and their mixture generalized CDTT model f(R) = R + + aR 2 , where 
a = ±1 a is a positive real number and fi is a parameter with units of mass. 
Both the CDTT and generalized CDTT models allow de-sitter solution with 
Ro = VSfJ 2 , whenever /i > 0. By choosing \i ~ H , the present accelerat- 
ing expansion can be achieved and even early time cosmological evolution 
can be recovered. Since laws of gravity gets modified on large distances in 
f{R) models, this allows several interesting observational signatures such as 
modification to the spectra of the galaxy clustering [3j H] , cosmic microwave 
background (5J E] and weak lensing [7J E\ ■ 

In order to provide f(R) gravity as a consistent theory of gravity, f(R) 
models are severely constrained cosmologically and by gravitational physics 
[9]- |TT] . The best known scale of gravitational physics is the solar system 
scale. To have stable stellar configuration, the function f(R) should satisfy 
the condition, d 2 f/dR 2 > at each point inside the star [12]. If this condition 
is not fulfilled, the star would collapse to some new configuration due to 
small perturbation. Also, this theory should compared to parameterized 
post-Newtonian (PPN) approximation, where the non-trivial metric elements 
outside the stars are g 00 = 1 — ^ and gn = 1 _ lppN %™ with observational 

value | 7pptv — 1 |«C 1. However, most of the j{R) models gives ^ppn = |, 
contradicting the observations. 

Gravitation theory and relativistic astrophysics have gone through exten- 
sive developments to the discovery of stellar collapse. A star has a life cycle 
wherein it is born in gigantic clouds of dust and galactic material, then evolve 
and shine for millions of years and eventually enter the phase of dissolution 
and extinction. It happens when internal pressures subside by nuclear pro- 
cess, gravity takes over, and the star begins to contract and collapses onto 
itself. Due to the achievements of f{R) gravity in numerous areas of cosmol- 
ogy and astrophysics, it is quite natural to discuss dynamics of gravitational 
collapse in this theory. In recent years, we have explored some aspects of 
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dark energy and gravitational collapse in f(R) theory [15]- [21]. 

The problem of dynamical instability is closely associated with formation 
and evolution of self-gravitating objects. The pioneer work in this direction 
was done by Chandrasekhar [22]. Afterwards, this issue has been investi- 
gated for adiabatic, non-adiabatic, anisotropic and shearing viscous fluids by 
Herrera et al. [25J-[2I]- m general, dynamical instability discussed in terms 
of adiabatic index T. The expression for isotropic spheres may be given by 
T > | + n—, where n is a number of order unity that depends on the struc- 
ture of a star, M and R are the mass and radius of a star respectively. For 
example, for white dwarfs n = 2.25. Instability limit of collapsing stars also 
depends upon the physical quantities describing the properties of different 
fluids filling the stars. For example, the dissipating quantities increases the 
instability range at Newtonian corrections but makes the fluid less unstable 
at relativistic corrections [28]. Also, Chan et al. found interesting results by 
studying the effects of anisotropy, radiation and shearing viscosity at New- 
tonian (N) and post-Newtonian (PN) regimes [29]- [31]. Sharif and Azam 
[321 133] . have discussed the effects of electromagnetic field on the dynamical 
instability of spherical and cylindrical symmetric gravitational collapse. We 
have investigated the effects of a well known f(R) = R + aR? model on the 
dynamical instability of expansion-free spherically symmetric gravitational 
collapse [33] . 

It is a fact that most of the collapse models analyzed so far are spherical 
due to their wide astrophysical significance. However, the study of non- 
spherical collapse remains a major uncharted territory, even many attempts 
have been made, e.g., [351 ESj- Specifically, final fate of collapse of a non- 
spherical cloud coming from numerical relativity and certain analytical solu- 
tions in cylindrical symmetry provide some new examples about gravitational 
collapse [57] 158]. Nevertheless, a comprehensive analytical treatment of col- 
lapse still needs to be done. In this connection, we would like to extend our 
work on dynamical instability to the cylindrically symmetric collapsing stars 
in f(R) gravity. 

The manuscript is laid out as follows. In Section 2, we provide a full de- 
scription of the matter distribution, the line element, both inside and outside 
the fluid boundary, and the field equations. We also formulate the dynamical 
equations which governs the dynamics of gravitational collapse. In section 
3, we consider CDTT model and discuss its feature on cosmology and grav- 
itation. Also, we apply perturbation scheme to all the function and f(R) 
model. Considering Newtonian and post-Newtonian limits instability range 
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in the form of T factor is obtained. In the last section 4, work is summarized 
followed by an Appendix. 



2 Dynamical Equations 

The modified form of EH action leading to f(R) gravity can be written as 



S = ~L] d ' x ^9f{R), (2.1) 

where k is the coupling constant and f(R) is any arbitrary function of Ricci 
scalar. The field equations in metric approach are reached by variating above 
action with respect to g a p as follows 



F{R)R a p--f{R)g a p-V a V p F{R)+g aP nF{R) = kT^, (a, (5 = 0,1,2,3), 

(2.2) 

where F(R) = df(R)/dR, □ = VV,, with representing the covariant 
derivative and is the standard minimally coupled stress-energy tensor. 
Alternatively, this can be written in a form similar to the field equations of 
General Relativity (GR) 

G a ji = p[TZp + TZp], (2.3) 

where 



rpc ^ 



f(R) RP(R) gap + Va V,F( J R) - 9ttfi DF(R) 



(2.4) 



This way of writing field equations help to study the dark side of the universe 
in term of curvature terms provided that it does not satisfy the usual energy 
conditions. Consequently, this theory may be used to explain the expansion 
of the universe and other aspects of dark energy in the gravitational physics. 

We consider a cylindrically symmetric collapsing star bounded by a hy- 
persurface S. The interior region inside the boundary can be represented by 
the following line element 

dsl = A 2 (t, r)dt 2 - B 2 (t, r)dr 2 - C 2 (t, r)d6 2 - dz 2 . (2.5) 
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In order to preserve cylindrical symmetry, following constraints may be im- 
posed on the coordinates 



oo < t < oo, 0<r<oo, 



oo < z < oo, < 9 < 2n. 



(2.6) 



The exterior region across the boundary of that star can be represented by a 
cylindrically symmetric manifold in the retarded time coordinate v as follows 
[39] 

2M(» 



ds 1 



-dv 2 + Idrdv - r 2 (d6 2 + ~f 2 dz 2 ), 



(2.7) 



where 7 2 



~, A being the cosmological constant while M(u) is the to- 



tal mass inside the boundary surface. As gravitational collapse is a highly 
dissipative process, we assume that fluid filling the collapsing cylinder is dis- 
sipating energy in the form of heat flux q and can be represented by the 
following energy-momentum tensor [281 E3] 



(2.8) 



T a/ 3 = (p + p)u a up - pg a/3 + q a U/3 + q$u c 



where p stands for energy density, p for pressure and u a for four-velocity of 
the fluid. In co-moving coordinates, these quantities satisfy the relations 



u 



A' 1 5%, u a u c 



1, q a = qB-^* 



U q D 



0. 



(2.9) 



The effects of dissipation describe a wide range of situations. For example, 
using quasi-static approximation, limiting cases of radiative transport have 
been studied in [101 . It is found that hydrostatic time scale is very small as 
compared to the stellar lifetimes for different phases of a star's life. It is of 
the order of 27 minutes for the sun, 4.5 seconds for a white dwarf and 10 -4 
seconds for a neutron star of one solar mass and 10 km radius |41j . 

For interior metric, Eqs. fl2.3p yield the following set of field equations 



B'C 
~BC 



C" 
~C 



ABC 
\a + B + C 



C A'C 
~C + ~AC 



BC 
~BC 




pA 2 + 



A 2 (f-RF 



f" f' a 

B 2 ~ + B 2 ~~C 



A' R 
- qAB + - [ F' - --F - —F' 



A 



B 



(2.10) 
(2.11) 
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F C F' 
+ MC + B 2 \ A 



A" B AB MB' 
AB \ ~B ~ A + ~M ~ ~W 



F" F B F' M 
B 2 + MB ~ ~FP~A 



(2.12) 



K 

F 



1 

P 

K 



f 



-RF F 
~2 



(2.13) 



A" 



B 



+ 



AB MB' 



AB 2 A 2 B MB AB 3 
MC BC k 
= " F 



AB 2 C A 2 BC 



1 

p + - 

K 



AC 
MC' 

f 



C 
MC 
RF 



B'C 



C" 



B 3 C B 2 C 
F _ F" 
+ M~B 2 



FIB C\ F' ( M C r 
+ M \B + CI~B 2 (^ + ~C 



(2.14) 



Here dot and prime denote derivatives with respect to t and r respectively. 
Subtracting Eq. fl2.13p from f |2. 1 2j) . we have 




B' 
~B 




C 
M 



(2.15) 



The dynamical equations help to investigate the evolution of gravitational 
collapse with time and yield the variation of total energy inside a collapsing 
body with respect to time and adjacent surfaces. We have formulated these 
equations by using contracted Bianchi identities both for the usual matter 
and effective energy-momentum tensor carrying curvature terms. These are 
given by 



(m) (c) 

rpa{3 _j_ rpa0 



( M (c) \ 

u a = 0, T^ + T^ Xq = 0. (2.16) 
The two independent components of Bianchi identities for the system under 
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consideration read 




+ M r >t) 



0. 



(2.17) 



(2.18) 



where 0i(r, t) and <f>2(r,t) are given in Eqs. fl5.ip and ( 15. 2 p respectively in 
Appendix. 



3 CDTT Model and Perturbation Scheme 



where a was taken equal to —1 in the original paper [2]. It is claimed that 
with the parameter \i chosen as the order of the inverse of the universe age, 
\T X ~ 10 18 sec ~ (lO^eV) -1 , this model can be applied to cosmology and 
describes the present accelerating expansion of the universe. However, this 
additional term was not essential to describe earlier time inflation. Further, 
CDTT model has a homogeneous solution of de-Sitter metric with constant 
scalar curvature R = Rq = \^3fi 2 . Thus, in order to have the desired physi- 
cally plausible late time behavior, we have to set a/3/U 2 ~ H$. 

Henttunen, et al. [42J have analyzed the properties of polytropic stars 
by considering CDTT model with the conclusion that the density profiles 
in general resemble the Newtonian Lane-Emden solutions. They analyzed 
that interior solution for metric components and curvature of such stars are 
always singular, consequently, f(R) = R — fi 4 /R model is not experimentally 
suitable to describe the space-time around the sun. A possible way to avoid 
is to relax the requirement data set for the central boundary conditions, but 
a more plausible approach is to modify the functional form of f(R)- As 
it is mentioned in the introduction that sign of f"(R) determine whether 
this modified theory approaches the GR limit or not at high curvature. For 
f"{R) > 0, the model behave very close to GR and is stable and vice versa. 
Due to this reason, the original CDTT model was severely criticized [151 WQ . 



In this section, we consider the following form of CDTT model, 



f(R)=R + a^, 



(3.1) 
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In order to avoid instability, Hu and swaicky introduced CDTT model 
with a = 1 [32]. Thus, throughout the paper, we shall use following form 

W) = R+£. (3-2) 

For the present model, the stability condition is satisfied with extremely small 
and positive f"(R) everywhere inside and near the polytrope. Also, it aims 
to explain the current accelerated expansion of the universe. At least for this 
model, Einstein evolution can be recovered for most of the cosmic history, 
especially standard epoch of matter domination can be obtained, providing 
a sufficiently long time to satisfy observations. It is mentioned here that the 
additional inverse term will not significantly alter standard evolution until 
today and that the solution lies well within present constraints from Big-Bang 
Nucleosynthesis [IB] . 

It is noticeable that CDTT model adds a perturbative function in the 
form of inverse power to the EH action. Thus for such a modification, the 
Einstein cosmology may be recovered in the limit that the perturbation term 
disappear. This system has great similarity with perturbation theory [47J. 
First, we analyze our work on dynamical instability when /i ^ 0, then we 
reduce our results asymptotically, fj, — > 0, being the Einstein solution. 

In this paper, the perturbation scheme introduced is used to analyze the 
instability conditions of the dynamical equations. We assume that initially 
all the function whether material or metric have only radial dependence, i.e., 
lie in the static equilibrium. However, afterwards, all these quantities have 
time dependence in their perturbation. Taking < e C 1, quantities may 
be written in the following manner 



A(t,r) 


= A {r)+eT{t)a{r), 


(3.3) 


B(t,r) 


= B (r) + eT{t)b{r), 


(3.4) 


C(t,r) 


= rB(t,r)[l + eT(t)c(r)}, 


(3.5) 


p(t,r) 


= Po(r) + ep(t,r), 


(3.6) 


p{t, r) 


= pa(r) + ep(t,r), 


(3.7) 


m(t, r) 


= mo(r) + erh(t, r), 


(3.8) 


q(t,r) 


= eg(t,r), 


(3.9) 


R(t,r) 


= R {r) + eT{t)e{r), 


(3.10) 


f(R) 


= [R (r) + 2/i 4 i? V)] + eT(t)e(r)[l - 2 f i i R^(r)}, 


(3.11) 


F(R) 


= 1 - /i? - 2 (r) + 2e^R-,\r)T(t)e{r). 


(3.12) 
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The static configuration of the field equations (j2.10j) - (l2.12p and ( I2.15P is 
obtained as 



-[Po + 



K 



{l + 2V(-3i?o^o +<)}] 



1 


[(I) 2 - 


^0 


*0 


Bl 




r5 


5 



(3.13) 



1 A' n (B' 



PrO 



fi 4 Rn 



K 



_ K 2fi i R ' i R [ 
A B 



4 D-3 d/ 


1 - ^Rf 



2R o 2 R'o 
B 2 



Bo r 



A B r 



2% 
rBo 



-DO 



(3.14) 



(3.15) 



Applying static part of the perturbed system, the first dynamical equation 
( 12.17P is identically satisfied while ( 12. 18[) yields 



A' 

Pr + (PO + Po)~T + K 4>2s = 0, 



(3.16) 



where 2s is the static part of equation (f> 2 . The perturbed configurations of 
Eq. fl2~T7D leads to 



,Ao 
Bo 



A 


f2A' 


Bo 


\A 


'2b 


c 

+ B 


Bo 



Bo 
Bo 



>ip 



0. 



(3.17) 



Here <pi p denote the perturbed form of equation <\>\ and is given in Appendix. 
Applying perturbation on Eq.f l2.lip and eliminating q, we have 



+ 



kA B 





1 

B~o 



— + b' + c' + ^(b + c) 
r A 



R7 4 



l-^Ro 



A' 



2^(e' + AeRo'R'o - e-£) f + — 



f. (3.18) 



Differentiating above equation and substituting in the second and third term 
of Eq. (I3.17P and then integrating resulting equation with respect to "t" , we 
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have 



P 



.(2b c 

ipO+PrO) I + 



(3.19) 



where 03 (r) is given in Appendix. Considering the second law of thermo- 
dynamics, we can express a relationship between p and p r as the ratio of 
specific heat by assuming an equation of state of Harrison- Wheeler type as 
follows 0H1 S9] 

p = T Pr0 p. (3.20) 

PO + PrO 

Here T measures the variation of pressure for a given variation of density. We 
take it constant throughout the region that we want to study. Substituting 
value of p from Eq.f l3.19p in f)3.20p . we have 



i> ' ( 7r • 77 ) ; 1 

. -DQ -DQ . 



Po 



Po+Po 

The perturbed configurations of second Bianchi identity is 



(3.21) 



B A' 

p +q^ + (p + P)^ + (po + Po) 



a 

T 



T 

T H 2p 

K 



(3.22) 



where <p2p represents the perturbed part of 02 and is provided in the Ap- 
pendix. Substituting values of q, p and p from Eqs.f l3.18j) . (I3.19P and (I3.2ip 
respectively in the above equation, we have 



+ 
+ 



T, 

K 



(2b c\ po 

-Po -5- + -5- + : <? 

\B B J po + Po 

T, N ( 2b C 

(Po + Po + I>o) "B" + -rT 

V-DO -DO 



T + (p + Po) 



a 



r- 



Po 



Po +Po 



+ 1 



T 
T 



>2p ~ 

Rq 4 



p 4 R ' 



kAq 



1 

~B~o 



— + b' + c' + ^(b + c) 
r A 



P 4 Ro 



-2 



2p\e' + 4eR^R' - e^] 

A o 



rBn 



f = 0. (3.23) 



It is mentioned here that above equation is the required evolution equation 



10 



for further analysis. The Ricci scalar curvature is given by 

A" MB' 




R = -'"T + AB 

B AB C 
B + AB + C 



MC B'C r 
~AC + ~BC 

AC BC 
~AC + ~BC 



(3.24) 



The static part of the above equation is obtained as follows 



-2 



A" 
An 



B n 



A B 



+ 



M n 



rAa rBn 



BL B' n 



Bl 



(3.25) 



The perturbed configuration of the Ricci scalar curvature with the use of 
Eq. fllEZgD yields 



aAZ 

An 



2bM ' 

~b7 



rBl 



A B 



b' + c'B' 



cB' 3bB' 



B 



aA' B' 



3bM B> 

Bn 



+ 



2c" 
rBl 



+ 



B^r 2 



b 

B~ 





( b 


2c\ 




)- T f 


\rB 


r J 


if 



This equation can also be written as 

f{t) - Mr)T(t) 



0. 



(3.26) 



(3.27) 



where 04 (r) is given in the Appendix. In order to find the instability range, 
we assume that all the terms of equation 04 are such that it remains positive. 
Consequently, the solution of Eq. (13.271) is obtained as 



T(t) = -e 



4>it 



(3.28) 



In the following subsections, we use above value in the dynamical equation 
Eq.f l3.23]) to analyze the evolution of the physical variables at the surface of 
the star. 



Newtonian Limit 

In this regime, we assume that p ^ Po an d A = 1, B = 1 in order to 
fulfill the Newtonian limit. Thus substituting M = 0, B' = in Eq. (13.23]) . 
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we obtain 



1 - /i 4 i?o 2 



2& + c + y + _ 2//i? 4 (e' - 4e.Ro ^ 1 



+ 



r{-p (26 + c)} 1 + a / po + 



02p(AT) 



ft' 



1 - ^ 4 i?Q 

r = o. 



(3.29) 



It is mentioned here that 4>2p(N) represent terms belonging to the Newtonian 
regime of perturbed second Bianchi identity. Substituting value of T from 
Eq.f l3T28|) . we have 



r < 



a'po + 05 0) 



M2& + C)]" (3 - 30) 

where 05 (r) is given in Appendix. It is noticeable here that adiabatic 
index depends upon the energy density, pressure and curvature terms in the 
Newtonian regime. Thus the collapsing system would remains unstable until 
Eq. CTIj) satisfied. 



Asymptotic Behavior 

When /i — > 0, the expression for T becomes 

i'po + ^(^ + b' + c') 



r< \po(2b + c)]' 

This results represents the GR solution. 



(3.31) 



Post Newtonian Limit 

Here, we include relativistic effects upto order m& . In this case, we have 

A = l-^, B = l + ^, (3.32) 



Ai 







m 



B' 



o 



-m 



Aq r(r — m ) ' B r(r + mo) 
Substituting above values in Eq. fl3.23p . we obtain 

^ [(p o +Po)(26 + c)-03(P^)] +(p +Po) [a(l + *?)]' + Mr 



(3.33) 



r > 



\p (2b + c)(l - ff) - + [p (2b + c)(l - - fftoai ' 



P0+P0 J 

(3.34) 
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In the above expression Mpn) denote those terms of 3 which belong to 
PN regime whereas 06 ( r ) is given in Appendix. Thus, system would be 
unstable in PN-approximation as long as the above inequality is satisfied. It 
can be observed that how relativistic and curvature terms are affecting the 
instability range of collapsing star. Further, in order to fulfill the dynamical 
instability condition, we need to keep all the terms positive. Hence, we 
assume that all the quantities involved in the above expression are positive. 
In addition, following constraints should be satisfied 



(PO + Po) (26 + c) > 3 (P7V) , 

p (2b 



-w, "0, P0<P3(PN) 

c)(l ) > — -. 

r Po + Po 



(3.35) 
(3.36) 



along with ^ < 1. 



Asymptotic Behavior 

As /i — > 0, adiabatic index T has the same representation. However, in this 
case 03 and 6 respectively reduced to 



M r ) 



Mr) = ~ 

K 



mo 

TK 



K 



2b + C I rj I — / I bm 

h o + c H — 



r 

2b + c 



1 f 2b + c . bm 

- + b' + C 1 + 



+ b' + t 



2 f2b + c 



+ b' + d 



3b + c 3m : 



+ ^ + (b' + c')[l + 



2m 



(3.37) 
(3.38) 



It is mentioned here that in both the N and PN regimes heat flux is not 
affecting limit of dynamical instability. 



4 Summary 

In this paper, we have assumed dissipative fluid configuration of a non-static 
cylindrically symmetric collapsing star and discussed the dynamical insta- 
bility for the final stages of the star by a perturbative approach. As study 
of gravitational collapse could be used as a paradigm to understand stel- 
lar formation. We introduced this issue in modified f(R) theory of gravity. 
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We have considered CDTT model which leads standard cosmic history. As 
short-timescale instabilities are avoided with extremely small and positive 
f"(R), so we have considered modified CDTT model with positive term 
It is worthwhile to mention here that this model satisfying the condition 
f"(R) > have stable high curvature limits and a well behaved cosmo logical 
solution with proper era of matter domination. 

To see the effects of CDTT model on the dynamical instability of fluid evo- 
lution during the collapsing process, we have applied a perturbation scheme 
on the field equation and dynamical equations. It is worth mentioned here 
that recently we have worked on spherically symmetric gravitational collapse 
evolving under the expansionfree condition [31]. In this study, the range of 
instability is independent of adiabatic index T showing the consistency of 
obtained results with the expansionfree condition which requires that fluid 
would evolve without compressibility. However, no such condition is imposed 
in this work, hence results appeared in the form of adiabatic index. 

Also, let us make the comparison of our results with Chandrasekhar's 
results on the limit of dynamical instability of spherically symmetric config- 
uration of matter [22] • In this work, dynamical instability depends upon the 
numerical value of the adiabatic index. It is found that if T > 4/3, the pres- 
sure in a star is strong enough than the weight of the outer layers which make 
a star stable. Whereas, for T < 4/3, the weight increases very fast than the 
pressure and star collapses resulting a dynamical instability. As concerned 
to work done in this paper, results depends on the physical quantities, like 
energy density, pressure, curvature terms and mass of the cylinder. Thus, in 
Newtonian regime collapsing star yield the instability limit in the the form 
of Eq. fl3.30p . By reversing the sign stability of the system would be achieved. 
However, in PN regime instability range not only depends on Eq. fl3.34p but 
also on certain constraints in Eq. fl3.35p and 03.36 p . 

Finally, it is mentioned here that our previous work had been done on su- 
per gravity model while present study made on CDTT model, hence the focus 
of future work is the combination of both the above models, i.e., generalized 
CDTT model in any of the symmetry. 



14 



5 Appendix 



0i M) 



+ 



+ 



+ 



X 



+ 



+ 



+ 



+ 



'>2s = 



K 



A 


f- 


kA 




B 


' F 


~^B 


A* 


C 


F 




A* 




ABC 
A + B + C 



+ 



f' a 

A^~C 



3A' B' a 
~ + ~B + ~C 



+ 



A 2 



K 



A 2 B 2 



F' - —F - —F' 
A B 



(5.1) 



<h{r,t) = 



B 2 

K 



1 



A 2 B 2 
f-RF 



F>-*F-*F> 
A B 



A' 
kA 

kB 

c_ 
1 



IB 2 
—F 



F F C F' 
+ A^ + A 2 C + B 2 



'A' B' V 
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